Storage of light: A useful concept? 
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We show both analytically and numerically that photons from a probe pulse are not stored in 
several recent experiments. Rather, they are absorbed to produce a two-photon excitation. More 
importantly, when an identical coupling pulse is re-injected into the medium, we show that the 
regenerated optical field has a pulse width that is very different from the original probe field, ft is 
therefore, not a faithful copy of the original probe pulse. 
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' Recent experiments [1,2] on the extremely slow propagation of a pair of optical pulses in a resonant medium and 
7-H , the subsequent coherent regeneration of the "probe" pulse have generated much excitement in the fields of optical 
computing and quantum optics. The physical explanation offered to these extraordinary phenomena rests on the 
adiabatic theory developed a decade ago by Eberly and coworkers [3]. To a large degree, the "dark-state polariton" 
theory [4] and all other recent works formulated around it are basically a field theoretical reformulation of the adiabatic 
approximation of pulse pair propagation in resonant media originally studied reasonably complete in the early 90's. 
These new treatments have provided some further understanding of the process, especially from the view point of 
^C) " atomic spin wave excitation and atomic coherence. Two main results of the these works, however, do not, in our view, 
accurately describe the physical picture of such a highly nonlinear excitation process. In this letter, we show, based 
on the adiabatic theory, that the claim of the "revival of the probe pulse" in recent experiments [1,2,5] is not accurate. 
Indeed, we show that under the conditions of these experiments the original probe photons are neither "stopped" or 
■ "stored". Rather, they are absorbed to produce two-photon excitation. Furthermore, we show that the " retrieval of 
<*S || the probe pulse" is not achieved without restrictions being placed on the characteristics of the second coupling pulse 
i 1 „ . used to regenerate the "probe" field. In fact, we show that the regenerated field has time-spectra characteristics that 
are very different from the original probe field. Therefore, "storage of light" is not an accurate description of the 
process. In addition, we show that all experimental observables can be well predicted with the usual treatment of 
combining classical electrodynamics and the three-state model without invoking a full field theoretical methodology. 
This is not surprising since when a detailed handling of spontaneous emission is not critical and the fields are not 
extremely weak, a mean-field approximation to the quantum treatment of the electromagnetic field is valid. Our 
analysis and results on the physical picture are unique and provide a clear understandings of the physical process. To 
the best of our knowledge the subtle understandings presented here on the "tapped light" problem are not contained 
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in either the original studies or any of the subsequent studies, including all recent theoretical works, on the subject. 

We start with a A system (Fig. 1) with the assumption that a probe laser (E p , frequency uj p ) is tuned on or near 
resonance with the |1 >^ |2 > transition. In addition, a coupling laser (E c , frequency u> c ) is tuned so that exact 
two-photon resonance between states |1 > and |3 > is achieved. We therefore have the following atomic equations of 
motion 
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dt r 



iVL p A 2 , (la) 



^ = in* p A! + in* c A 3 + % (s + C^) a 2 , (ib) 
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= i£l c A 2 , (lc) 



where Aj and jj are the jth amplitude of the atomic wave function and decay rate, respectively, and we have 
made phase transformation to remove z dependent phase factors. Also, &^,t c \ — ^21(23) ^p(c)/(2?i), S = uj p — io 2 \, 
lu p — uj c = k>3i, and t r = t — z/c is the retarded time which, in vacuum, is the combination of t and z that the 
laser field amplitudes depend on. Notice that we have included non-vanishing one-photon detuning < 5 < il c , a 
feature that is not included in both the original adiabatic theory and those developed recently around the "dark-state 
polariton" theory. 



In order to correctly predict the propagation of the probe and coupling laser pulses in a resonant medium, Eq.(l) 
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must be solved simultaneously with Maxwell equations describing the propagation of the probe and coupling fields. 
For plane waves in the slowly varying amplitude approximation the Maxwell equations resume the form 



dz 



in 12 A\A 2 , (2a) 



dfl* " 

~di) t = "*«-^- (2,,) 



where ^12(32) = ^NcOp^D^^l 2 1 '(fie). The process involves first seeking, from Eq.(l), an adiabatic solution for 
the atomic response, and applying the result to Eq.(2) to self-consistently describe the propagation of the pulse pair 
within the adiabatic limit. Following this procedure, we thus obtain adiabatic solution to the atomic response [3] 

' 4 ><-''-» - wxf- (3a) 

it \ i dAi i dA 3 

-5* <3b) 
- -WM- (3c) 



where fl(z, t r ) — y/\£l p (z, t r )\ 2 + \£l c (z, t r )\ 2 . As usual, with the adiabatic approximation used in a resonance situation, 
we must have fi c already strong when Cl p starts to build up. One scenario where the adiabatic approximation would 
hold through the whole laser pulse is if the two lasers peak at the same time, but the pulse length of the coupling 
laser is much longer. This feature of the pulse lengths and the requirement \ fl c r\ » 1 through the pulse length, r, 
of the probe laser will make results based on the adiabatic approximation quite accurate. 



With Eq.(3) as atomic response, we now solve Eq.(2) that now resumes the form of 

' dn Uc)\ _ Ki2 ( 32)T d ( rn; {c) 



dz J tr v/lrOcI 2 + \tQ p \ 2 dt r /T \^\tQ c [ 2 + \tQ p \ 
We first note that the quantity 



(4) 



F(Mr) = M + M (5) 

represents the sum of the photon fluxes at ui p and uj c divided by the concentration of the medium through which the 
waves propagate. Differentiate F with respect to z while holding t r fixed and apply Eq.(4), we immediately reach the 
conclusion that F depends only on t r . This permits one to evaluate F by evaluating it at z = 0, the entrance to the 
atomic vapor cell, e.g. 

f( , = o, t)= mr + mr, (6) 
K32 «i2 

Therefore, whenever F(z, t r ) occurs it can be replaced by F(0,t r ), as determined in Eq.(6)[6]. The lack of dependence 
of this quantity on z when the full set of equations is solved numerically is an important test for the validity of the 
adiabatic approximation. 

Let us consider, for simplicity, the case where K\i = n 32 - Since ^/|£! c | 2 + |£! p | 2 does not depend on z when t r is 
held fixed, Eq.(4) can be recast into 

K 12(32)T d ( %(c) T \ 




\n cT \ 2 + \n pT \ 2 ou/t \ y/\n cT \ 2 + \n P T\- 

Define 

w p = /T7 ^z= ==,w c = j= ;; c ; = = , m 




V(t r ) = / [\n c (0,t r ) T \ 2 + \n p (0,t r )T\ 2 d ^ },U(Z)= / K 12 Tdz , (8b) 



3 



Equation (7) now becomes 



where general "travelling wave" type solutions are immediately obtained as 

W p = F p (v-u), (10a) 
W c = F c {v-u). (10b) 

The functions F p and F c are easily determined by evaluating at z — (so u = 0). When the second coupling laser 
pulse is injected into the medium after a time delay, the predictions about the revival of the "probe" pulse are all 
contained in this solution. The tabulation of F p remains the same as long as only a second coupling pulse is sent into 
the medium at a later time. Notice that if |£! p (z, i r )| 2 << |$l c (z, i r )| 2 , the population in state |3 > is always small 
and the coupling laser propagates as in vacuum. In this limit, with the replacement of £l c (z, t r ) = fi c (0, t r ), Eq. (10) 
can still be used to determine Q p (z,t r ) even for Ku ^ «32- This is the essence of the approximations applied in the 
works by Rcf.[4]. 

We now examine the wave propagation and "probe revival" with the field profiles [7] 

n p (0,i) = n p0 e- ( r )2 , (lla) 

n c (o, t) = n M ( e -°- 2 (T) 2 + R e -°^-*o)^ (llb) 

Here, fl p o and Sl c o are real constants characterizing the peak amplitudes of the two half-Rabi frequencies before the 
pulses enter the resonant medium, R is the ratio of the Rabi frequency at which the coupling laser recurs to its initial 
amplitude, and x = ta/r is the value of t r /r at which the peak of the coupling laser recurs. 

We first consider the case where R = and the group velocity of the probe pulse is sufficiently small so that the 
coupling laser dies away before the probe pulse can propagate through the cell. When the coupling laser begins to 
die out, the intensity of the two lasers becomes proportional to each other, therefore the relation 



v{t r ) ~ \n c0T \ VW2 + \n p0 T\ 2 ^/^/2 (12) 

is appropriate. Also, |O p (0,i r )| 2 has long been very small compared with |f2 c (0, t r )\ 2 . Thus, for t r /r > 3.0 we have, 
as a very good approximation 



n* p {z,t r ) ~ r!(0,i r )F p ((|O c0 r| 2 yW2+ \fl p0 T\ 2 y^j2) - k 12 tz). (13) 

At such a late time during the pulse, the argument of F p depends only on the z coordinate so that the time dependence 
of CI* is obviously exactly the same as that of Cl(0,t r ). When |fi c (0, t r )| is several times larger than |f2 p (0,t r )|, as 
always is at such late times at z = 0, this means that Cl*(z, t r ) has the same time dependence as fi*(0, t r ) at such late 
times. Note that at late times Cl*(z,t r )/Cl(z,t r ) = As(z,t r ) is also independent of retarded time, indicating that as 
Q p (z,t r ) approaches zero, the ratio of populations in |1 > and |3 > stays fixed. It is precisely this way of having two 
fields to go to zero with fixed ratio that preserves the the adiabatic approximation during the process. There must 
be a coherent superposition of states in |1 > and |3 >, such that A^{z,t r )/A\{z,t r ) = Cl*(z,t r )/fl*(z,t r ). A similar 
relation holds when the second coupling laser pulse starts to build up. In order for the behavior to remain adiabatic, 
the probe laser half-Rabi frequency must build up proportional to the coupling laser, with the ratio of the two being 
the local ratio of A3/A1. This persists until the depletion of the population of |3 > forces the ratio to decrease at 
later t r . By choosing z such that 2n\2TZ = |O c ot| 2 a/57t /2 + \J ir / 2\Q p vt\ 2 , we thus have the same argument of F p at 
z = and t r = t = 0, resulting F p = il p o/ y/\fl c o\ 2 + |^ P o| 2 - This is the largest value that F p takes on, and at this 
depth into the medium the value of 1^3(2, t r )\ matches its largest value at z = 0. This population persists at large t r 
until very slow collisional effects destroy the coherence left behind in states |1 > and |3 >. It is this long persistence 
of a coherent mix of populations in states |1 > and |3 > that leads to the regeneration of an optical field when a 
delayed second coupling pulse is injected into the medium. Of course, the number of atoms left in |3 > is (within the 
adiabatic approximation) equal to the number of photons in the original probe pulse. 

We now investigate the case where £I c qt = 20, fl p oT — 5, 72T = 0, K12T = K32T = 200cm -1 , and R = 4. Based 
on our adiabatic theory, we predict a maximum population of \A$\ = 1/Vl7 at a depth of 2.86 cm into the medium. 
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This is indicated by the line of constant color in a contour plot of |-A3(z,t r )| (see Fig. 2a) leading from t r = and 
z = out to the horizontal path at z ~ 2.86 cm. A corresponding surface plot for f2 p r as functions of t r /r at 
z = 3 cm is given in Fig. 3a which shows the long asymmetric tail on f2 p as described above. In the region between 
2.5 < t r /r < 5 the ratio of the two half-Rabi frequencies is close to constant, averaging around 0.24. This ratio is also 
close to the adiabatic approximation for A 3 since |£! p | 2 << |£l c | 2 . In the same region, however, there is no optical 
field left, indicating no photons are "stored" or "stopped". Every probe photon is converted to the excitation of the 
state 1 3 > [8]. 

We now come to an important point, e.g. to show that the regenerated field has characteristics that are very 
different from that of the original probe pulse. That is, the regenerated field is not a faithful copy of the original 
probe pulse unless careful restrictions are placed on characteristics of the second coupling laser pulse. Recall that 
F p was determined from the functional dependence of the probe and coupling laser pulses at z = 0, in particular, 
W p (v(t)) = — Az(0,t), we thus have, during the recurring coupling laser pulse 

v(U) =S+ ^ c0 r| 2 yW2 (l + erf ^^2/5 { ±-M.^ , (14) 

where S = \Q c ot\ 2 -^/5tt/2 + \Q p0 t\ 2 Equation (14) is of central importance in the following analysis on the 

characteristics of the regenerated field. We first note that R must be large enough so that v(t r ) — Ki 2 rz m > 0. That is, 
R must be large enough so that the group velocity of the regenerated photons are large enough to exit the cell before 
the laser induced transparency ends. Three time markers, therefore, are important for describing the regenerated 
field when it reaches the end of cell where z = z m . The first marker is the earliest time at which F p (0) = 0. This 
time marker is determined by v{t r \) — n\2Tz m = 0. At a later time the value of F p (S/2) will be equal to —^(0,0). 
This second marker represents the time when F p , for the given set of parameters, reaches it maximum value at z m , 
and it is determined by v(t rm ) — K\<2jz m = 5/2. Finally, the third marker is the time at which the regenerated pulse 
completes its exit from the cell, i.e. the time at which F p (S) — and v(t r2 ) — n\2Tz m = S. When these relations are 
used in Eq.(14), we immediately obtain 

n;( Zm ,t rm ) = aw-^-**)'/^) m fe= , 2 - (is) 

Notice that if R is chosen to make the argument of F p at the time t r = td exactly the same as its value at z — and 
t r = t = 0, we then have 

\n p0 T\ 2 ^7/8 + \n c0T \ 2 v/W8(i + R 2 ) = n 12 Tz m , (16) 

and 

n*(z m ,t d ) = flifl c oi . p0T (17) 

By determining the value of t r such that v(t r ) — Ki2TZ m = or S, a range of time over which the regenerated pulse 
rises from zero and returns to zero at the exit of the cell may be obtained. We thus estimate the FWHM pulse length, 
in the unit of the original probe pulse length r, to be 



_VW2/' 1 



tt p o 



a 



cQ 



(18) 



The key results shown in Eqs. (14-18) indicate that the regenerated field is not a replica of the original probe pulse. 
Therefore, the concept of "storage of light", in the reported experimental studies [1,2], is not accurate. Furthermore, 
Eq.(15) clearly indicates that the instantaneous phase of the regenerated field is closely related to that of the second 
coupling pulse, in additional to the phase of the atomic coherence P32, therefore, cannot be just that of the original 
probe pulse alone. In fact, if the probe laser has only penetrated a small fraction of the thickness of the vapor cell 
when the coupling laser pulse has passed by, then R will turn out to be much larger than unity. This means that the 
width of the regenerated field is generally much smaller than the width of the initial probe pulse. Correspondingly, 
the bandwidth of the former will be much larger as can be seen from Fig. 3a. Both Fig. 2a and 3a show that in 
"reviving the probe pulse" every regenerated photon comes at the expense of flipping population from state |3 > to 
|1 >. When the population of |3 > has been exhausted, there can be no further photon generated. 
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Extensive numerical calculations carried out by simultaneously solving Eqs.(l-2) have shown very good agreement 
with the adiabatic calculation described above. In Figs. 2b and 3b we show a contour plot of | A 3 (z, t r )\ and a surface 
plot of Q* for the same parameters given in Figs. 2a and 3a. From Fig. 3b, we notice that there is no probe field 
between t r /r = 2 and the arrival time of the second coupling pulse, i.e. there is no photon "left" or "stored" in the 
medium. All probe photons have been absorbed in producing the coherent excitation of the state |3 > as can be seen 
from Fig. 2b. The case with non-vanishing one-photon detuning in our adiabatic theory also produces a result that 
is in very good agreement with numerical calculations. Numerical calculations have also been vigorously tested by 
making use of the fact that with 7 2 r = 0, the sum of the squares of the three state amplitudes should be unity. In 
all numerical examples described in this paper the condition of unity was preserved through at least seven significant 
figures if 7 2 t = 0. 

We have shown analytically and numerically that when < S < |O c |, a coherent optical field with the frequency 
very close to that of the original probe can be regenerated by re-injecting a coupling pulse. Detailed analysis on the 
conditions and characteristics of this regenerated field, including the estimate of the pulse width, has shown that it 
is not the replica of the original probe pulse. We, therefore, caution the use of the concept of "storage of light" in 
the context of recently reported experiments, since there has no evidence that this is indeed the case [9]. In fact, 
one can couple the states |1 > and |3 > with a coherent magnetic pulse to create the coherence required. Under 
this circumstance, if an optical pulse that couples the states |2 > and |3 >, commonly referred to as the coupling 
pulse, is injected into the system, an optical pulse that couples the states |2 > and |1 >, commonly referred to as the 
probe pulse, will be generated. The latter is most certainly not the replica of the magnetic pulse used to create the 
coherence. 
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6. Equation (7) can also be derived exactly from Eq.(3) for the case of slowly varying |A 2 | 2 and not too large 72. 

7. These field profiles are reasonable description of the experiments reported. In fact, as long as the adiabatic 
conditions are observed different profiles will lead to the similar conclusions. In addition, numerical simulations 
have suggested that even in the region where the adiabatic conditions are in valid, a similar conclusion can still 
be reached, provided certain non-adiabatic restriction are observed. 

8. Figure 3a shows that all probe photons are converted to the excitation of the state |3 >. If one can regenerate, 
from the excitation of the state |3 >, photons with optical properties identical to that of the original probe filed, 
then it would be correct to claim the revival of the probe field. As we will show later, however, this is not the 
case in recent reported experiments. 

9. We emphasize that the purpose of the present study is to show that photons from the probe pulse are neither 
"stopped" or "stored" in the recently reported experiments [1,2]. Our intension is to stimulate further studies 
for truly storage of the quantum information. 



Figure captions 

Figure 1. Energy level diagram showing relevant laser excitations. 

Figure 2. Contour plot of A3 (z, t). a. Adiabatic solution, b. Full numerical solution. Parameters used: f2 p r = 5, 
fl c T = 20, 72T = 0, ki 2 t = k 32 t = 200cm~ 1 , R = 4, td/r = 11. 



Figure 3. Surface plot of Q p (z,t). a. Adiabatic solution, b. Full numerical solution. Parameters used: fl 
c t = 20, 7 2 r = 0, k 12 t = k 32 t = 200cm~ 1 , R = 4, t d /r = 11. 
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